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Abstract 

Elastic wave characteristics of the heat conduction in low-temperature thin wires can be studied 
via a wave scattering formalism. A reaction matrix formulation of heat conductance modeled by 
elastic wave scattering is advocated. This formulation allows us to treat thin wires with arbitrary 
surface disorder. It is found that the correlation in the surface disorder may significantly affect the 
temperature dependence of the heat conductance. 
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I. INTRODUCTION 



In a thin wire where the electron or phonon wavelength is in the same order of its width, 
a continuous description of transport is no longer valid. Instead, a quantized unit of electric 
charge e or a characteristic scale of heat energy ksT {ks is Boltzmann constant and T is the 
temperature) is needed to understand the transport properties. When such a wire is under an 
external field, e.g., a voltage difference in the electron case or a temperature difference in the 
phonon case, system-independent features such as the universal electric conductance quanta 
Ge = 2e^ /h {h is the Planck constant) and the universal heat conductance G = Air'^klT/Sh 
quanta will emerge With the universal heat conductance quanta successfully demon- 
strated recently 2l], there have been great interests in the wave nature of heat transport at 
the nanoscale. For example, recent experiments showed that the thermoelectric properties 
of silicon nanowires are much better than that in bulk silicon ^. 

Besides a nanoscale geometry, one needs low temperatures (e.g., T ~ 1 K) to observe 
a phase-coherent transport and hence the wave nature of phonon transport {3, 0, S]- At 
low temperatures only acoustic phonon modes are populated, with their characteristic 
wavelength much longer than typical atom-atom distances. As such, an elastic wave ap- 
proach becomes appealing for a quasi-classical description of heat transport. Indeed, elastic 
wave propagation in thin wires has been analyzed in several studies with fruitful results 



0,0, Big, 



10| . However, in these earlier publications only very simple geometries of a thin 



wire were considered, leaving the case of a thin wire with surface disorder unexplored. This 
motivates us to develop a framework that can handle elastic wave scattering for arbitrary 
wi.e geometHe. Specifically, we apply a reaction matnx tonnulat.on, previously thoroughly 
developed for electron transport in nanoscale waveguides [11], to the case of heat transport. 
We believe that this is the first time that a reaction matrix formulation is applied to the con- 
text of elastic wave scattering as well as heat conductance. In so doing we use the so-called 
thin-plate approximation, which means that one of the dimensions is decoupled from the 
other two. This thin-plate approximation is widely used in analyzing transport behaviors 



121. 



Our formalism here for treating elastic wave transport can handle arbitrary surface ge- 
ometry or arbitrary surface disorder. This feature can be very useful because, depending 
on the process of growth, there are often different types of surface disorder introduced to 
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nanowires 



13|. Preliminary studies on the effect of surface disorder on nanowire heat con 



16|- 



duction show the existence of universal features in the presence of disorder 
Nevertheless, it is still an important problem to understand in detail how surface disorder 
affects the actual temperature dependence of heat conduction. Furthermore, as learned 
from our previous work on electron transport in rough waveguides, conductance properties 
of a nanoscale rough waveguide can depend strongly on the involved energy scale and the 
long-range correlation of the surface disorder. We hence expect, similar to early studies of 
heat transport in one-dimensional systems with partially random defects or partially random 



coupling constants 



17, 
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19| . that long-range correlation in the surface disorder of thin 



wires may have important implications for heat transport properties. 

In addition to presenting details of a reaction matrix formulation of elastic wave scattering 
in two-dimensional geometries, we also discuss several specific results. For example, it is 
found that the throughput can be indeed modified by assuming different kinds of correlations 
in the surface disorder of a wire. In the one-mode regime the curve of transmission versus 
phonon frequency may develop a clear dip, thus considerably affecting the temperature 
dependence of the heat conduction. 

This work is organized as follows. In Sec. II we introduce our elastic wave scattering 
model for heat transfer in long thin wires. In Sec. Ill, we first introduce how specific surface 
roughness can be generated, and then show the details of how the throughput of out-of- 
plane elastic waves can be calculated using a reaction matrix formulation. Representative 
numerical results are presented and discussed. Section IV concludes this work. 



II. HEAT TRANSPORT MODELED AS A PROBLEM OF TWO-DIMENSIONAL 
ELASTIC WAVE SCATTERING 

Heat transport in a solid is a result of unbalanced phonon population. There are mainly 
two different perspectives in understanding the heat current in nano- materials [2^. One 
is based on the Kubo formalism, where phonon current is understood as a response to an 
external temperature field, with the proportionality coefficient expressed as a heat conduc- 
tivity tensor. Similar to the electron transport case, the conductivity tensor is related to 
the current-current correlation function. Because the current-current correlation function 
should be calculated for different phonons and because the phonon number can change (um- 
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FIG. 1: (a) Schematic plot of a out-of-plane mode of an elastic wave propagating in a thin plate, 
(b) Schematic plot of a two-dimensional rough waveguide that models a thin wire with a rough 
surface at y > 0. An, Bn, Cn and Dn denote scattering amplitudes. See the text for details. 

klapp process), the Kubo formalism leads to nonlinear equations that are difficult to solve 
in practice. The other perspective is due to Landauer, where the current can be determined 
at the surfaces of the sample by applying the correct boundary conditions. In particular, 
the known solution for outside a sample can be matched at the boundary using a scattering 
formalism. In this latter framework, from the right bath with temperature Tr to the left 
bath with temperature T^, the heat energy flux J is given by 



where w is the phonon frequency, Ga is the throughput of the mode a. Here tir^l is the 
phonon distribution given by rji{w) = l/[exp{hw/kBTi) — 1], where Tj are the temperatures, 
with indices i = R,L. Heat conductance is defined by energy flux per temperature difference, 
i.e.. 



where for small temperature difference (Tr — Tl), k, is approximated by the temperature 
derivative of the energy flux. After substituting rii{w) into Eq. ([1]) we have. 




(1) 



K = J/ AT ^ dJ/dT, 



(2) 




(3) 



4 



where we have defined the average temperature T = (T^ + Tr)/2. In a perfect wire through- 
put is unity for each open channel, then k in such a perfect case is given by 

^-^N, (4) 

where N is the total number of open modes. As seen from Eq. ([3]), the throughput Ga{w) is 
the main quantity that determines the temperature dependence of k in Landauer's formalism. 
In the language of a scattering problem, Ga{w) is determined by the absolute value squared 
of certain scattering amplitudes. 

To calculate the throughput of elastic waves in a thin wire, we use a thin plate approx- 
imation to treat full elastic wave equations. A schematic plot of a two-dimensional rough 
waveguide modeling a nanowire with surface disorder is shown in Fig. 1(b). Because at 
low temperatures contribution from any optical modes is exponentially small we use an 
elastic theory to describe elastic wave scattering in a relatively long wire. As we show in 
Fig. 1(b) as an example, we intend to solve a scattering problem in a wire with the width- 
length ratio set to be 1:100. We scale all the lengths by the wire width, which is taken to 
be unity. 

The elastic wave equation for our model system is given by 



dt 
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ciV'u+{cf-ci)V{V-u), (5) 



where u is the displacement vector , q is the longitudinal phonon velocity and q is the 
transverse phonon velocity [3]. For convenience we also scale the velocities by q. That is, 
in our calculation q = 1. In the thin-plate approximation, the gradient in the 2;-direction is 
assumed to be zero, thus simplifying Eq. ([5]). Such an approximation decouples the in-plane 
acoustic modes from the out-of-plane acoustic modes. Assuming that the in-plane and out- 
of-plane modes behave in a similar fashion, we confine ourselves to the out-of-plane modes 
only [see Fig. 1(a)]. Technically speaking, implementing the reaction matrix formulation 
for out-of-plane modes is also simpler than for in-plane modes. The wave equation for 
out-of-plane modes becomes 

1 d'^u,. d^u^, d^u^ 

+ (6) 



cf dt"^ dx"^ dy"^ 

where Uz{x) is the displacement in the z direction. The boundary condition for this wave 
equation will be specified below. 
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The next step is to calculate the scattering of the elastic waves in a two-dimensional 
geometry as shown in Fig. 1(b). We will employ a reaction matrix formulation for this 
scattering problem. The basic idea of a reaction matrix approach is as follows. First, the 
system is divided into the lead region and the scattering region. Second, the known solution 
in the lead region is matched at the boundary with the found solution in the scattering region. 
This matching is implemented by using a basis set satisfying the appropriate boundary 
condition. In the following section we will show how the reaction matrix can be constructed 
based on appropriate basis states. 



III. SCATTERING THROUGHPUT IN A THIN WIRE WITH SURFACE DIS- 
ORDER 



As learned from electron scattering in a rough waveguide [22|, it can be expected that 
disorder- induced localization will also play an important role in elastic wave scattering in a 
thin wire. Furthermore, if the localization length is smaller than the wire length, then the 
scattering throughput Ga{w) is exponentially small; if the localization length is much larger 
than the wire length, then Ga{w) should be close to unity. Because in the low-temperature 
regime only a few scattering channels contribute to heat transport, it is also important to 



realize that the localization length might depend strongly on the phonon frequency |16 |. 
In particular, based on the Born approximation, the electron localization length is inversely 
proportional to the structure factor (defined below) of the correlation function of the surface 
disorder evaluated at twice of the scattering electron wavevector [22, l23|. By direct analogy 
between the Schrodinger equation for electron matter waves and the elastic wave of matter, 
we can expect that the correlation function of thin-wire surface disorder can also affect 
strongly the localization length of a scattering elastic wave. Therefore, we are interested 
in (i) if it is possible to manipulate the frequency dependence of the scattering throughput 
of elastic waves, and (ii) to what extent different surface disorders might affect the heat 
conductance. 

For completeness we first discuss how to computationally produce a random surface 
characterized hj y = 1 + ri{x) [see Fig. 1(b)]. This can be done by dividing a long wire into 
pieces, shifting each piece up or down randomly, and then connecting each piece smoothly 



with a cubic spline function 22, l2J]. Randomness of the surface function r] thus generated 
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can be described by an auto-correlation function {rj{x)ri{x')) = cr^Cnix — x'), where a is 
the variance of r]{x). The Fourier transform of this auto correlation function is defined as 
the structure factor, denoted xi^)- For disorder close to a white-noise type or for disorder 
without long-range correction (computationally, we can only divide the wire region into a 
finite number of pieces, e.g., 100 pieces, hence the surface disorder we generate is at most 
close to a white-noise type), the structure factor trivially takes a constant value over a broad 
range of k. To generate a more realistic rough thin wire with "colored" or long-range surface 
disorder, we define another function p{x), 



sin ( a*" 


X 


)- 


sin( 


x) 




a'' 


X 



where and a'^ are two introduced parameters. Consider then a new function r]{x) resulting 
Tom a convolution of ri{x) with p(x), i.e., fj{x) = J dx' p{x')r]{x — x') . It can be easily shown 



22] that due to this convolution, x(^) iio"^ becomes a "square bump" function, with the 
bump edges located at a' and a^. That is, x{^) is zero for k < a! oi k > a!' and is a constant 
in between. This convolution method is used below to produce several "bump" functions 
as an "engineered" surface-disorder structure factor, with different bump widths. Certainly, 
one may also combine a number of different p{x) functions to generate a rather arbitrary 
structure factor x(A;). 

A. Constructing the reaction matrix 

The reaction matrix method is a well-known time-independent scattering formalism. As 
such we need to consider the stationary solutions to the wave equation in Eq. (El). Let us 



first substitute the ansatz exp(— ztyt)^^ into Eq. ([6]) 2l|], yielding 



=0 (8) 



2l| 



where w is the phonon frequency. The realistic boundary condition for is given by 

8^ = "- 
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A stationary solution to Eq. ([8]) for a straight channel satisfying the above "zero-stress" 
boundary condition at ?/ = and y = 1 can then be found, namely, 

^ni = - i?n^ cos(^); (10) 



V 

where and are the scattering amplitudes. The wavevector is given by 



kn= \/w^ — {ji'K I dy, where n = 0, 1, 2, ■ ■ ■ . Note that when fc„ takes imaginary values, the 
solutions are called "evanescent modes". As mentioned above, in our unit system d = 1, 
the wire length is L = 100, and the phonon frequency w is in terms of the value of the first 
mode, i.e., Wi = Ctii/d = vr. We also scale temperature by T^, where Td = Tiwi/ks- For 
a thin wire of a width d = 10~^cm and for an acoustic transverse phonon speed (Silicon) 
Ct = 5.8410^ cm/sec, we have To = 14K. 

In a reaction matrix formulation of wave scattering, solutions in Eqs. (11 01) and ffTTl) can 
be taken as "free" solutions in the left and right leads that are connected by the wire. To 
find the relation between the wave amplitudes An, Bn, Cn, and Dn, these two solutions will 
be matched with that in the wire. Hence, we need to solve the elastic wave equation in Eq. 
(IH]), confined by two boundaries defined hj y = l + fj{x) (or y = l+ri{x) for an almost white- 
noise type surface disorder) and y = [see Fig. 1(b)]. To that end a set of eigenf unctions of 
Eq. (IE]) satisfying the zero-stress boundary condition are needed. In addition, as one main 



feature of the reaction matrix formulation 
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22l |. the derivative of the eigenfunctions with 



respect to the scattering coordinate x should be zero at the lead-wire interface, i.e., at x = 
and X = L. 

To find the eigenfunctions in the scattering region, we employ a coordinate-transformation 
technique. That is, we consider two new coordinates x' = x and y' = y/[l + fi{x)], where 
fj{x) is the surface function of the wire. In the x' — y' coordinate system, the elastic wave 



equation 



Refs. 
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becomes complicated. However, it is now considered over a straight channel (see 



22l | for details of this coordinate-transformation technique). In the new coordinate 



system an eigenstate 0„ satisfying the necessary boundary condition can be written as 

n^\y') = $^/5;,,^,^=i_cos(^)cos(mV), (12) 

m,m' V \ ' / V / 

where m, m' are integers, the summation over m and m' must be truncated for practical 
purposes, and are the expansion coefficients to be found numerically. 
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Once the eigenstates with eigenfrequencies Wi are obtained, we impose the con- 

tinuity condition at a; = and x = L, yielding a reaction matrix, 

R....tmm, (13) 

^-^ W — Wi 

i=l 

where M is the size of the involved basis sets (exact for M = oo), and 0^ = (jf cos{mTy) 
is the overlap of ith eigenstate with the rith "channel" solution cos(n7ry). The index n,m 
of the reaction matrix Rnm denote the scattering channels associated with real wavevectors 
kn and km- The reaction matrix Rnm can then yield the scattering matrix, 

I-iKR 

* = TTTkr' 

where / is a unit matrix of the same size with R, and K is a. diagonal matrix with diagonal 
elements related to the wavevectors kn 



24l |. The scattering matrix then relates the 




incoming wave with the outgoing wave via 

Finally, the total scattering throughput is defined as G = X]q — J2nm\^n,m\'^, where 
the summation is over all open modes. G essentially describes the total transmission of a 
scattering wave up to a certain given phonon frequency, without taking into account the 
actual thermal distribution of these phonon modes. 



B. Throughput and heat-conductance results 

At low temperatures, very few acoustic modes can be populated. Hence it is of interest 
to examine the scattering throughput in the one-mode regime. This is shown in Fig. 2. 
Specifically, Fig. 2(a) depicts the results associated with two realizations of a wire with 
almost white-noise surface disorder. It is seen that when the phonon frequency is small, 
the throughput G is fluctuating drastically, with its average value far from unity. This 
indicates that in these two cases the localization length of the elastic wave is comparable 
to the system size. As the phonon frequency increases, G is seen to be close to unity, 
suggesting that the localization length is now well beyond the wire length. Figure 2(b) 
shows what happens if we impose a correlation on the surface disorder, but with the noise 

9 




FIG. 2: Acoustic phonon transmission in the one-mode regime, as reflected in the dimensionless 
throughput G. (a) Results for two reahzations of a waveguide with almost white-noise surface 
disorder, (b) Results for three different cases of colored surface disorder, with the associated 
structure factor set as a "square bump" function with three different bump widths. See the text 
for details. 

variance fixed at a = 0.08. In particular, the dotted curve (curve 2), the dashed curve 
(curve 3), and the solid curve (curve 1) are for "engineered" surface functions generated 
from a convolution of a random surface function r]{x) with [sin(0.77rx) — sin(0.47ra;)]/ (O.Tttx), 
[sin(0.77rx) — sin(0.57rx)]/(0.77rx), and [sin(0.67ra;) — sin(0.57rx)]/(0.67rx), respectively. It is 
seen that G develops a clear window for certain phonon frequencies. Despite statistical 
fluctuations, these throughput windows also show differences in their widths, namely, curve 
2 shows the widest window and curve 1 shows the narrowest window. This is consistent 
with the differences in the bump widths of their associated surface structure factor. We 
also examined other realizations of random surfaces with the same surface structure factor, 
confirming that the location of the throughput windows shown in Fig. 2 does not change 
much. We have also checked that the location of the throughput window quantitatively 
matches the profile of x(2fc), the surface structure factor evaluated at twice of the wavevector. 

In Fig. 3 we examine how different properties of the wire disorder might affect the heat 
conductance k calculated from Eq. (3). To make a better connection with the results 
in Fig. 2 we examine a temperature regime where thermal excitation only populates at 
most two channels. First of all, for the case without disorder (upper thin solid curve), 
the plotted quantity k/T is seen to be roughly a constant for low temperatures. This is 
somewhat expected because a perfect thin wire is known to show a linear temperature 
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0.15 0.2 0.25 0.3 

FIG. 3: The ratio between heat conductance k and temperature, as a function of temperature, for 
wire geometries considered in Fig. 2. Thin soHd hne is for a wire without surface disorder. Line (4) 
is for a wire with almost white-noise surface disorder. Other hues are for wires with their surface 
function generated from a convolution of white-noise with some smooth function explicitly given 
in the text. Here and in all other figures below, the plotted quantity is always dimensionless and 
Ko = kswi. 

dependence in k at low temperatures. Certainly, as the temperature increases, deviations 
from the linear behavior emerge due to an increasing contribution from the second populated 
channel. By contrast, in the presence of surface disorder, the linear temperature dependence 
of K might not hold from the very beginning of the temperature range shown in Fig. 3 (curves 
2,3,4). Interestingly, for different surface disorder, k/T shows quite different behavior. The 
almost white-noise case (curve 4) deviates most from the case without disorder, with the 
maximal deviation around a factor of two. These observations suggest that one might 
achieve some subtle control over the temperature dependence of k in the low-temperature 
regime by manipulating surface correlations, or say something about the surface correlation 
by measuring k at different temperatures. 

To shed more light on how the throughput windows observed in Fig. 2(b) affect the 
temperature dependence of k, let us now consider a model throughput curve that allows for 
analytical calculations. Consider first two step functions /a(w) and fhiw), namely, faiw) = 
ior w < a and fa{w) — 1 for w > a; and fb{w) — for w < b and fb{w) — 1 for w > b. 
We denote the heat conductance as M(a) for G{w) — fa{w) and M(6) for G{w) — /&(w), 
respectively. Assuming that there is only one open mode, we obtain from Eq. (3) that 

H(a) = ksw, (^p^T) - 21og (^"^^ - 1) a + ^ - 2r5ft [Li, (e«/^)]) (16) 
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FIG. 4: (a) Heat conductance from direct numerical integration (dotted) of Eq. (3) versus the 
analytical result (solid line) of Eq. (jl6p . for a = 0.2. As indicated by the inset, the dotted line is 
on top of the solid line, (b) Same as in (a) but the dashed line is calculated from an expression in 
Ref. [l]. Note the different scales between panel (a) and panel (b). 

and an analogous expression for K(6), where is the so-called dilog function. Interestingly, 
at this point we found that the widely cited Ref. l|] does not contain the second term in Eq. 
f fT6|) and misses a factor of two in the third term. To ensure that our expression is indeed 
correct we show in Fig. 4(a) that direct numerical integration results based on Eq. ([3]) are 
identical with the analytical result of Eq. ( |T6l) for a = 0.2. For the sake of comparison we 
also show in Fig. 4(b) the associated analytical result of Ref. jl| (dashed line), as compared 
with the same solid line in Fig. 4(a). Note also that because the function Li2{z) for certain 
values of z can be explicitly evaluated, it is now possible to predict the exact value of K(a) 
for certain e"'^'^ . 

Let us now examine the heat conductance behavior if G{uj) is a square-well function 
Gsw{w)^ i.e., it is zero for a < w < 6 and is unity elsewhere. Because Gsw{w) = fbiw) + 1 — 
fa{w), one obtains from Eq. (fT6|) 



K 



+ K(6) - K(a) 



where T is already scaled by T^. Approximately, for b > a » T, we have 

T\b-a) 



K = h ksWi 



{T-a){T-b) 



and for a < 6 << T we have 



K 



+ ksWi [ a — b + 2a log(a) — 2b log(6) 



b' 



36r2 



(17) 



(18) 



(19) 
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FIG. 5: Temperature dependence of k/T, for throughput curves modeled by square-well functions. 
The arrow indicates an increasing width of the square-well function, in steps of 0.08, ranging from 
0.0 to 0.4. Note that quantities here are plotted in the logarithmic scale. 

These two limiting results show that in either case [Eqs. f|T8l) or Eq. f|T9|) ]. a window in 
the throughput curve can cause k/T to decrease as compared with what is expected from a 
linear temperature dependence of k. 

From Eq. (fT7|) it is also possible to examine how the width (6 — a) of the square well 
function Gsw{w) as a throughput curve affects k/T. To that end we show in Fig. 5 a 
logarithmic plot of k/T versus T/Td. The arrow therein indicates an increasing {b — d) in 
steps of 0.08, which ranges from 0.0 to 0.4. Interestingly, as the throughput window increases 
its size, k/T decreases. The k/T function is also seen to deviate from a constant function 
more and more, and a better description of the temperature dependence of k seems to be a 
power function. We note that this observation is consistent with experimental observations 
in superconductive materials {4]. Note also that the analytical result here also explains the 
general trend seen in our numerical results in Fig. 3. That is, the lowering of k/T curves in 
Fig. 3 is associated with the widening of the throughput windows shown in Fig. 2(b). 

So far we have studied the low-temperature regime where one or two modes are signif- 
icantly populated. What happens if more phonon modes are populated? To answer this 
question we show in Fig. 6 the throughput Giw) for up to five modes, for a wire without 
disorder [case (a)], with almost white- noise surface disorder [case (b)], and with colored sur- 
face disorder [cases (c) and (d)]. A number of observations are in order. First, G[w) can be 
nonzero for phonon frequencies much smaller than Wi, a feature different from that in elec- 
tron transport. In the latter case one can only excite the first mode after going beyond the 
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FIG. 6: Throughput G{'w) versus phonon frequency. Case (a) is for a perfect wire, case (b) is a 
wire with white-noise surface disorder, and cases (c) and (d) are for two reahzations of a wire with 
colored surface disorder. For a better view cases (c), (b), (a) are shifted upwards by two, four, and 
six units, respectively. The arrows indicate the dip in the throughput curves. For example, the 
right arrow for curve (d) indicate a throughput dip in the two-mode regime. 

associated cut-off energy. This difference arises from the fact that the Schrodinger equation 
for electrons and the elastic wave equation for phonons have different boundary conditions. 
Second, in the presence of surface disorder, the threshold frequency for opening a new chan- 



nel and hence a drastic increase in G 
This is analogous to the electron case 



w) gets higher as compared with the noiseless case. 



22| |. which can be roughly explained via a decrease in 



the effective wire width. Third, for cases (c) and (d), the dip in G{w) is clear only for the 
one-mode and two-mode regimes. For higher modes, the throughput dip is no longer clear 
and is buried by fluctuations in G{w). 

Can those throughput curves in Fig. 6 manifest themselves differently in the heat con- 
ductance behavior, when the temperature is high enough to populate many phonon modes? 
The answer is positive based on the results shown in Fig. 7. In particular. Fig. 7 depicts 
the temperature-dependence of k, for those cases shown in Fig. 6. It is seen that these cases 
still show quite different behavior even when five phonon modes are significantly populated. 
Cases (c) and (d) associated with colored surface disorder lie between the curve for a perfect 
wire and that for a wire with almost white-noise surface disorder. Results in Fig. 7 indicate 
that for temperatures of the order of ~ 10 K, different surface disorder properties of a thin 
wire of a width ~ 10 nm may still cause a considerable difference in heat conductance. The 
results are also plotted on a log-log scale in the inset of Fig. 7. However, we do not see 
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FIG. 7: Temperature-dependence of heat conductance for those cases considered in Fig. 6. Case 
(e) is for a perfect (straight) wire with twice the width, in which the temperature dependence of k 
is approaching the law for bulk materials. Inset shows the same data in a log-log plot. 

any evident power-law dependence of k versus T. This is somewhat expected because the 
T-dependence of k is quite complicated [see, for example, Eq. (17)]. For the sake of com- 
parison, a perfect wire with twice the width is also shown in Fig. 7 (e). In that case 
the temperature dependence is much stronger. This is anticipated because of the scaling 
of K in bulk materials. 

IV. CONCLUDING REMARKS 

In this work we adopted an elastic wave picture of heat transport at sufficiently low 
temperatures. We presented, for the first time, a reaction matrix treatment of the scatter- 
ing of elastic waves in thin wires with surface disorder, with realistic zero-stress boundary 
conditions implemented. We have shown that correlations in the surface disorder can lead 
to structures in the frequency-dependence of throughput and can hence affect considerably 
the temperature dependence of heat conductance. Though not done in this study, the heat 
conductance for a thin wire with an experiment-related surface correlation function can be 
examined in a straightforward manner. Because our treatment is a natural extension of early 
electron transport studies, insights gained from studies of electron waveguide transport can 
now be much relevant to understanding heat transport in thin wires. 

Because the boundary condition for in-plane modes of elastic waves is much more diffi- 
cult to handle for a thin wire of arbitrary shape, here we restrict ourselves to out-of-plane 
modes only. However, preliminary results show that in-plane modes have similar qualitative 

15 



behavior as out-of-plane modes. 
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